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Abstract 

The covariant Dirac equation in Robert son- Walker space-time is studied under the 
comoving coordinates. The exact forms of the spatial factor of wave function are 
respectively acquired in closed, spatially flat, and open universes. 
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1 Introduction 



The Robertson- Walker metric plays a central role in modern cosmology. The expanding 
universe is usually described by it. Therefore the behavior of relativistic particles 
obeying the covariant Dirac equation in Robertson- Walker space-time has evoked many 
investigations in history. Such works go back to Schrodinger [1] , Brill and Wheeler [2] , 
Barut and Duru [3]. A general discussion of the case of massless neutrino has been 
given by Brill and Wheeler. 

In literature, Li and Guo ever considered the exact solutions of Klein-Gordon equa- 
tion and the Dirac equation in Beltrami-de-Sitter coordinates [4]. They gave right 
solutions for Klein-Gordon equation and acquired the discrete spectrum [5] of the mass 
of a spin-0 particle because of the continuous condition of the Legend function at the 
boundary of de Sitter space-time. But they obtained the wrong evolutionary equation 
from the Dirac equation, thus gave the wrong solutions for the Dirac equation in de 
Sitter space-time. 

Had been aware of the importance of the behavior of a spin-| particle in the expand- 
ing universe, Barut and Duru investigated the exact solutions of the Dirac equation for 
three typical models of expanding universe [3] . But they assumed that the space-time 
is spatially flat. Thus the exact solution of the Dirac equation in general Robertson- 
Walker space-time haven't been given in literature. 

In this paper, we investigate the behavior of a massive spin-| particle in the open, 
spatially flat, and closed Robertson- Walker space-time, respectively. The exact forms 
of the spatial factor are obtained. Our work will be a correction and extension to the 
paper of Li and Guo, on the other hand it is also a complementary to the work of 
Barut and Duru. The exact solutions of the Dirac equation in Robertson- Walker will 
give a stage for studying the possible usage of neutrino oscillations [6] in cosmology, 
we think. 



2 



2 Dirac equation in Robert son- Walker space-time 

The cosmological principle, the hypothesis that all positions in the universe are essen- 
tially equivalent, assumes that the "smeared" universe is homogeneous and isotropic 
about every point. The convenient implement to describe this "smeared" universe is the 
Robertson- Walker metric, which is manifestly spherically symmetric. In the comoving 
coordinates, the Robertson- Walker metric of the universe is taken to be* 

dr 2 



ds 2 = g^dx^dx v = dt 2 - R 2 (t) 



+ r 2 (d6 2 + sin 2 6d V 2 ) 



(1) 



1 — kr 2 

where R(t) is an unknown function of time, which is called cosmic scale factor, and k 
is a constant, which by a suitable choice of units for r can be chosen to have the value 
+1, 0, or —1. When k = +1, 0, or —1, the space-time is respectively closed, spatially 
flat or open. 

To set the stage for solving the Dirac equation in a curved space-time, we introduce 
the local tetrad field as follows* 

<T = e^(x)e\(x)rj ab , (2) 

where g^ u being the space-time metric, and the metric tensor of Minkowski space-time 
Vab is 

r/ 00 = +1 , r] 11 = rj 22 = r] 33 = -1 , rf = for a ^ b . (3) 
In a curved space-time, the covariant Dirac equation reads 



* r(x) d^~ Zyi(x)rV(;r) 



tp(x) = m-0(x) , (4) 



Here 7 M (x) are the curvature-dependent Dirac matrices and r^(a;) are the spin connec- 
tions to be determined. The curvature-dependent Dirac matrices 7 M (x) are presented 
in terms of the tetrad field as 

Y(x) = e» a (xh a , (5) 



^we adopt the conventions that x° = t, x 1 = r, x 2 = 9, x 3 = ip and % = c = 1. 
■''In this paper, using Roman suffixes to refer to the bases of local Minkowski frame; using Greek 
suffixes to refer to curvilinear coordinates of space-time. 



where 7" denotes the standard flat-space Dirac matrices, which satisfies {7 s , 7^} = 2rf h . 
The tetrad field in the case of Robertson- Walker metric (1) is thus taken to be 

1 \ 



e» a ( x ) = diag ^1. 
Inserting Eq.(6) into Eq.(5) yields 

7°(x) = 70 , 
7 2 (x) = 



Vl - kr 2 1 



R(t) ' R(ty R(t)r sin 0, 



7 W = ^ — 7i 



(6) 



-72 , l 3 (x) 



R(t) 
1 



i?(t)r ' ' ' w i?(t)r sin 
The spin connections r^(x) in Eq.(4) satisfy the equation 



73 



(7) 



(8) 



where IT are the Christoffel symbols for the metric (1), which are determined by 



1 ur ( dgrp , <9g„ r 9e 



pp <->& I rv .. ~r 



pp 



dx^- dx p dx T I 
Hence Eq.(8) can be solved for the spin connections r^(x) which we determine as 

R(t) 



(9) 



r 



, r\ 
1 



2y/l - kr 2 



7o7i , 



1 



T 2 = - J R(t)r7o7 2 + -VI - kr 2 ^2li , 

I . 1 \ 

T 3 = -R(t)r sin 6»7o7 3 + -VI - kr 2 sin 6*737! + - cos 6*7372 , 



(10) 



where R(t) denotes — ^p, so that the combination 7 M (x)T M (x) in Eq.(4) simplifies to 



, . , . 3-R Vl — kr 2 cos6 
l^)Tp(x) = --70 + -^-7! + 2/?(t)rsin ^7 2 • 



Inserting Eq.(7) and Eq.(ll) into Eq.(4), then we obtain the Dirac equation in 
Robert son- Walker space-time 



1 d_ 

r sin 9^ 3 dtp 



r ^ l\ 1 (d COS 9" 



mi? } ip{t, r, 9, if) = . 



(12) 
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3 The exact solutions of the Dirac equation 

We can simplify the Dirac equation (12) by writing [4] 

ip(t, r, 9, (p) = (sin 9)-*R-*V(t, r, 9, <p) . (13) 
The reduced equation in terms of r, 9, <p) is of the form 



R + im^j 



v / T^^7o7i (^ + 1) + ^(6, tp) * , (14) 



where 



^'^ )= 707l ( 7 4 +73 iin^^) (15) 



is a Hermitian operator, as we know, which is defined in Ref. [7] first. The eigen equation 
of K is 

K^!^ = q% fi , ? = 0,±1,±2,..., (16) 

where the eigen value q is an integer, which has been proven in Ref. [2, 7]. 

Only the two matrices j and 7x remain explicitly in the simplified Dirac equation 
(14), they can therefore be represented by 2 x 2 matrices 

70 = ( I -°i ) ' 71 = ( ° T ) ' (17) 

where 1, Ti are respectively 2x2 unit matrix and the Pauli matrix, namely 

*-(;;)■ m) 

We separate the angular factor from the wave function, and represent the radial and 
temporal factor by a two-component spinor, that is 

*« = <MM ($;§). < 19 > 

where the angular factor, f (0, y?), is determined by the requirement K"0 f (0, y>) = 
<^O f (#, </?). This equation for eigenstates of the angular motion has been investigated 
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by Schrodinger [7] . He found, the operator K is related to the total angular momentum. 
Their eigenvalues satisfy [8] 



S = T(j + 



-(I + 1) for j — I + - 



I 



for j — I 

J 2 



(20) 



Also the eigenfunction 6 f is related with the spherical harmonics as follows 

Q<(0,<p) oc Y^je,^ oc P r ; i+ hcos6)e^ + ^ , 



(21) 



or 



e<(0,tp) oc Y^ifrip) oc Prh^6)e^-> , (22) 

where P(cos6>) being the associated Legendre functions, and m± \ — 0, ±1, ±2, • • •, ±Z. 

Using equations (16), (17), and (19), we can rewrite the equation (14) into two 
equations 



R \Wt + im ) 



<f>(r,t) 



d 



Vl-kr 2 (!- + -) +- rix(r,t), 



(23) 
(24) 



To solve the above equations, we separate the functions 4>(t,r) and x(^ r ) hito radial 
and temporal factors, respectively 



<j>(t,r) = C/i (r)Tx(t) , X (t,r) = U 2 (r)T 2 (t). 



(25) 



With the help of Eq.(18), inserting the above equations into Eq.(23) and Eq.(24) yields 
the evolution equations 



R {jt + im ) Tl{t) ~ i£T2 ^ = ° ' 
R \Jt ~ I ~ ieTl ^ = ' 



and the radial equations 



\ar r J r 
d 1\ c 



VI - kr 2 — + - 

\dr r ) r 



\ar r J r 
d 1 



Vl - kr 2 — + - +- 

\ dr r r 



(26) 
(27) 

U x + £ 2 ?7i = , (28) 
U 2 + e 2 /7 2 = . (29) 



Obviously, if U\ is substituted by U2, and simultaneously q are replaced by —q, then 
the equation (28) becomes the equation (29). Therefore, in the following, it is of only 
necessity to study the equation (28). 

The equation (28) directly simplifies to 



^(PUi (2 OI \dU 1 fqVl-kr 2 -q 2 , , w . 

The constant k appears in the above equation, the value of which will affect the exact 
solution of Eq.(30). We, therefore, investigate Eq.(30) according to the value of k 
respectively. 

3.1 In spatially flat Robertson- Walker geometry 

In the case of spatially flat Robertson- Walker geometry, the constant k satisfies k — 0. 
The equation (30) is thus reduced to 

ar z r dr \ r z J 

Through a simple transformation, one can easily find that the equation (31) is the 
Bessel's equation of order — q + \. The solutions of the above equation are [9] 

Ui(r) = ^=J ±{ ^ 2) {er) . (32) 



l Er - J -V> 2 

Where J ± ^_^{er ) are called the Bessel functions. 

3.2 In closed Robertson- Walker geometry 

In the case of k — +1, the Robertson- Walker space-time is closed. The equation (30) 
becomes 

^d 2 U x [2 \dU 1 (q^/T^-q 2 2 \ TT 
Do the transformation 



then, in terms of p, the radial equation (33) is rewritten as follows 



,d 2 f/i 1 



d\J x 



p( 1 -pY-nr + ^-p)(3-p)-^ + 



dp- 



dp 



2p 



Ap 



U ± = . 
(35) 



The above equation can be transformed into a hypergeometric equation, In the interval 
of < r < 1, the solutions of which can thus be represented as [9] 



U 1 (r,q>0) = 



vr 



i-i 



^1^72 + 1 



/ l l Vi-^-i 

• F r £+ 2'^ + 2 Vl-rHl 



C/i(r,<j<0) 



vr 



1+e 



/ i 3 vr^-r 



(36) 



Where F is the hypergeometric function. 

3.3 In open Robertson- Walker geometry 

In this subsection, we investigate the case of k = —1, where the space-time is open. 
The radial equation (30) is simplified to 
d 2 U, 



2 ,d z U 1 /2 



dUi 



dr \ r 
Do the similar transformation as Eq.(34), set 

VI + r 2 - 1 



We then rewrite the equation (37) as 



v / TT7 2 "+ 1 ' 



P(l-P) : 



(^ 2 )(l-p) ; 



2p 



4p 



+ e 2 + l 



(37) 



(38) 



C/i = . 
(39) 



Similar to the equation (35), the above equation can also be transformed into a hyper- 
geometric equation [9], the solutions of which are taken to be 



£/i(r,<r>0) 



l+ie 



£Mr,<r<0) 





(40) 



We have indicated that the solution of U 2 can be obtained from equations (32), 
(36), and (40) by substituting — q for <;. Hence we have acquired the exact forms of the 
spatial factor of wave function in Robertson- Walker geometry. 

In a different coordinates of spatially flat Robertson- Walker space-time, the authors 
of Ref. [3] have investigated the exact solutions of the Dirac equation for three models 
of expanding universes. It is obvious that our evolutionary equations (26) and (27) are 
equivalent to the equation (10) in Ref. [3]. Hence the solutions of equations (26) and 
(27) must be the same as the results given in Ref. [3]. 

The authors of Ref. [3] have given the exact solutions of three special models for 
expansion: 1. R — a t, a being a constant, a model of a linear expanding universe; 2. 
R = aoVi, a model of a radiation-dominated universe; 3. R = e Ht , H being a constant 
also, a model of de Sitter universe. These models are typical in cosmology, we will not 
discuss any more. 



spatially flat, and open Robertson- Walker space-time are given, respectively. This 
work is a complementary to the Ref. [3]. 

Acknowledgement: I would like to thank Prof. D. Du, Prof. C. J. Zhu, Prof. Z. 
Chang, Prof. A. Sugamoto and Prof. I. Oda for their help and encouragement. 



4 Conclusion 



The behavior of the spin-| particles in Robertson- Walker space-time is investigated. 
The exact expressions for the spatial factor of the spin-| wave functions in closed, 
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